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1. Introduction 

We are concerned with semileptonic decays of B mesons (B and B*) into orbitally excited P 
wave D mesons (collectively denoted as D**'s): B^*'^ D**lv. These decays are of particular 
interest, because there is a persistent conflict between theory and experiment, the so-called "1/2 
versus 3/2 puzzle": while experimental results indicate that a decay into "1/2 P wave D**'s" is 
more likely, theory favors the decay into "3/2 P wave D**'s" (for recent reviews cf. [|l|, ^). 

1.1 Heavy-light mesons 

A heavy-light meson is made from a heavy quark (b, c) and a light quark (u, d), i.e. 
B = {bu, bd} andD = {cu, cd}. 

In the static limit {mh,mc — > °°) there are no interactions involving the static quark spin. There- 
fore, it is appropriate to classify states according to parity and the total angular momentum of 
the light quarks and gluons j (cf. the left column of Table |l|). 

If mh,mc are finite, j is not a good quantum number anymore. States have to be classified ac- 
cording to parity ^ and total angular momentum J (cf. the right column of Table [l]). Although j is 
not a "true quantum number" anymore, it is still an approximate quantum number justifying the no- 
tation Z)j. The above mentioned P wave D**'s are {D* , D\ , Di , D*} = {D^^ , dJ^^ , Dp , 0^^]. 



f 


7^ 


(1/2)- ^ S 


= B,D 

1 = B\D* 


(1/2)+ ^ P 
(3/2)+ ^ P+ 


0+ = D* = Dq/^ 

1+ = d; = d|/2 

1+ = Di = dJ/^ 
2+ = D* = Dj/^ 



Table 1: Classification of heavy-light mesons (left: static limit; right: finite heavy quark masses). 

1.2 The 1/2 versus 3/2 puzzle 

Experiments (ALEPH, BaBar, BELLE, CDF, DELPHI, D0), which have studied the semilep- 
tonic decay B —>^ X^lv (where Xc is some hadronic part containing a c quark), find the following 
composition of X^: 

• « 75% D and D*, i.e. S wave states (which is in agreement with theory). 

• 10% Dj and D2 , i.e. j = 3/2 P wave states (which is in agreement with theory). 

1/2 

• For the remaining w 15% the situation is rather vague: a natural candidate would be Dq 

1 /2 

and Dj , i.e. j = \/2 P wave states. This, however, would imply 

1/2 3 /2 

r(B ^ 1 ' v) > T{B Zv), which is in conflict with theory. This conflict between 

experiment and theory is called the 1/2 versus 3/2 puzzle. 
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On the theory side most statements are made in the static limit mh,m^. 00. In this limit the 
eight matrix elements relevant for decays B —>■ D** 1 v can be parameterized by two form factors, 
the Isgur-Wise functions T1/2 and T3/2 1^- Here we only list two of these matrix elements: 



(D;/'(v')|c-757m^|B(v)) - Ti/2(w)(v-v')^ 

{DfiV,e)\cr,r,b\B{v)) oc T3/2(w)((w+i)£;;„v«-£;^vVv;), 



(1.1) 

(1.2) 



where v and v' are the four velocities associated with the B and the D meson respectively, 
w = (v' • v) and e is the polarization tensor of the D meson. 

By means of operator product expansion (OPE) a couple of sum rules has been derived in the 
static limit [Q, §]. The most prominent in this context is the Uraltsev sum rule, 



Jn) 
S/2 



(1) 



-Ski; 



1 



(1.3) 



where T1/2 = l'{/2' '^3/2 = '^3/2 ^^'^ '^^^^ all 1/2 and 3/2 P wave states respectively. From 

experience with sum rules one expects approximate saturation from the ground states, i.e. 



^3/2' 



1) 



h/2 



(1) 



1 

4' 



(1.4) 



which implies |Ti/2(1)| < I '^3/2(1) I - This in turn strongly suggests 



1 /2 J Z 

r(S — > Dq J Z v) < r(B Dj 2^ / v), which, as already mentioned, is in conflict with experiment. 

Phenomenological models ||^, 0] give the same qualitative picture, even when considering 
finite heavy quark masses [^]. 

Possible explanations to resolve the 1/2 versus 3/2 puzzle include the following: 

• The experimental signal for the remaining 15% of Xc is rather vague; therefore, only a small 

1/2 1/2 

part might actually be Dq and D/ . 



,3/2 



Sum rules like (1.3) might not be saturated by the ground states. 



Sum rules derived by OPE hold in the static limit and might change for finite heavy quark 
masses. 



• Sum rules make statements about the zero recoil situation (w = 1), where the B and the D 
meson have the same velocity; to obtain decay rates, however, one has to integrate over w. 

With a dynamical lattice computation of Ti/2(1) and 13/2(1) in the static limit, which is pre- 
sented in the following section, we attempt to shed some light on this puzzle. 



2. Lattice computation of Tj /2 and T3 /2 

For a more detailed presentation of this computation we refer to [Q]. We use a method, which 
was proposed and tested in the quenched case in [p^. 
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Since the "Isgur-Wise relations" ( pTT| ) and (|L^) are not directly useful to compute Ti/2(1) and 
'^3/2(1) (the right hand sides vanish at zero recoil), they have to be rewritten as shown in [[Til]: 



{D'J\v)\cy,yjD,b\B{v)) = -igju(m{D'J^) -m{B)y,i2{\] 
{Dy\v,e)\cYsyjDkb\B{v)) = +/V3£y,(m(D^/') 13/2(1) 
We compute T1/2 by means of ( ^?T| ) and an "effective form factor": 

'^^1/2(1) = lim Ti effective (^0 - -^2) 

'^l/2,effective(^0 — — ^2) = 



1 



(2.1) 
(2.2) 

(2.3) 
(2.4) 



To this end we need static-light meson creation operators ^t"^', iff'^^-^ and static-light meson 

masses m{S), m{PJ) and m{P^), 2-point and 3-point functions, and norms N{S), N{P-) and N{P+). 
Zcg is a perturbatively computed renormalization constant, whose derivation is explained in detail 



in [ |12| , g]. The computation of T3/2 is analogous. Explicit formulae can be found in [|9p. 
2.1 Simulation setup 

We use X r = 24-' x 48 gauge configurations produced by the European Twisted Mass Col- 
laboration (ETMC). The gauge action is tree-level Symanzik improved and the fermionic action 
Nf = 2 Wilson twisted mass at maximal twist yielding automatic ^(a) improvement of physical 
quantities. The lattice spacing is a = 0.0855 fm. To be able to extrapolate our results to physi- 
cal light quark masses, we consider three different bare quark masses /iq corresponding to "pion 
masses" mps, which are listed in Table ^ For more details regarding these gauge configuration we 



refer to 1 13, 14|. 





mps in MeV 


number of gauge configurations 


0.0040 


314(2) 


1400 


0.0064 


391(1) 


1450 


0.0085 


448(1) 


1350 



Table 2: Bare quark masses, pion masses and number of gauge configurations. 



2.2 Static-light meson creation operators 

The meson creation operators we use are latticized versions of the continuum expression 

^(r)(x) = (2(x)yt/nr(n)[/(x;x + rn)i//-W(x + m), (2.5) 

where 2(x) creates a static antiquark at position x, i//(")(x + rn) creates a light quark separated by a 
distance r from the static antiquark, ?7 is a gauge covariant parallel transporter and Y a combination 
of spherical harmonics and 7 matrices yielding well defined parity ^ and total angular momentum 
of the light degrees of freedom j. The operators are collected in Table ^. 
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r(n) 




f 


Oh 


lattice i'^ 


notation 


75 





(1/2)- 


Ai 


(1/2)- , (7/2)- , ... 


5 


1 


0+ 


(1/2)+ 




(1/2)+ , (7/2)+ , ... 


P 


71/11 - 72/12 (cyclic) 


2+ 


(3/2)+ 


E 


(3/2)+ , (5/2)+ , ... 


P+ 


75 (71 «i -72^2) (cyclic) 


2- 


(3/2)- 




(3/2)- , (5/2)- , ... 





Table 3: J: total angular momentum; j: total angular momentum of the light degrees of freedom; parity. 



2.3 2-point functions, static-light meson masses, norms of meson states 



With meson creation operators ( |2.5| ) at hand it is straightforward to compute the 2-point func- 



tions 



r G {75, 1, 7i«i -72«2}- 



(2.6) 



From these 2-point functions we extract the meson masses m{S), m{P^) and m{P^) via effec- 
tive mass plateaus. To illustrate the quality of our data we show effective masses for jUq = 0.0040 
in Figure |l]. For details regarding the computation of the low lying static-light meson spectrum 



within our twisted mass setup we refer to [[151, |16p . 



H^ = 0.0040 - m„„„i,,(S) and m,„,„,i,,(P_) 
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Figure 1: Effective masses for 5, P and P+ for /iq — 0.0040. 



Moreover, we obtain the ground state norms N{S), N{P-) and N{P^) by fitting exponentials 



to the 2-point functions (2.6) at large temporal separations. 

2.4 3-point functions 

The computation of the 3-point functions is again straightforward. We chose to represent the 
covariant derivative inside the heavy-heavy current in a symmetric way by a single spatial link in 
positive and negative direction. 

2.5 Results 



In Figure ga we show the effective form factors Ti/2,effective (cqn. (2.4)) and T3/2^effective for 



to — h = 10 as functions of to — 1[ for /iq = 0.0040 (plots for the other two quark masses look 
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qualitatively identical). We extract T1/2 and T3/2 by fitting constants to the central three data points 
as indicated by the dashed lines. Results are collected in Table Q 



0.6 

0.5 
0.4 
0.3 
0.2 
0.1 



= 0.0040 - T|;2,eltective 13,2 elfective (to - = 1 0) 



t 



I «■ 



„. ;;;;»;;i;i^irt-v::vts.v... 



i • I 



13,2 = 0.521 ±0.013 . 
T„2 = 0.300 ±0.014 . 



8 10 
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0.6 



extrapolation to the physical u/d quark mass (Iq - 12 = 10) 
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0.15 
- in GeV^ 



Figure 2: a) Effective form factors Ti/2.effective T3/2. effective for tQ — t2 = 10 and /Xq = 0.0040. b) Linear 
extrapolation of T1/2 and T3/2 in (mps)^ to the physical u/d quark mass. 





^1/2(1) 






0.0040 
0.0064 
0.0085 


0.300(14) 
0.313(10) 
0.309(12) 


0.521(13) 
0.540(13) 
0.524(8) 


0.181(16) 
0.194(13) 
0.178(9) 



Table 4: T1/2 ™d T3/2 and their contribution to the Urlatsev sum rule. 
As expected from sum rules T3/2 is significantly larger than Ti/2- Moreover, we find that the 



ground states fulfill the Uraltsev sum rule ( |1.3| ) by around 80%. 

We use our results at three different values of the pion mass to linearly extrapolate /2 and 
T3/2 in (?Mps)^ to the physical u/d quark mass (mps = 135MeV; cf. Figure ^5). Our final result is 



"phys 



"1/2 



(1) = 0.297(26) 



"phys 



"3/2 



(1) = 0.528(23). 



(2.7) 



3. Conclusions 



Our result (2/7) confirms the sum rule expectation that 13/2(1) '^1/2(1) in the static limit. 
When comparing to the experimentally measured form factors {T'^^j^{\) = 1.28 and '^^j^i^) = 0.75 



[]17[]) we find fair agreement for T3/2 but a strong discrepancy for Ti/2- 

In our opinion this discrepancy calls for action both on the theoretical and the experimental 
side: it would be highly desirable to have a first principles lattice computation of T1/2 and T3/2 be- 
yond the zero recoil situation and also for finite heavy quark masses; on the other hand a thoroughly 
refined experimental analysis of the decay into 1/2 D** 's, for which the signal is rather faint, seems 
to be necessary. 
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